% 事实验证2 欧拉-伯努利横梁
% 未完。。。


function A = Euler_Bernoulli(n)

%%%%%%%%%% get A
A = zeros(n);
A(1,1:4) = [16 -9 8/3 -1/4];
A(2,1:4) = [-4 6 -4 1];
A(n-1,n-3:n) = ([16 -60 72 -28]./17);
A(n,n-3:n) = [-12 96 -156 72]./17;

for i=3:n-2
    A(i,i-2:i+2) = [1 -4 6 -4 1]
end
%%%%%%%%%%